NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


19990921  070 

THESIS 


ANALYSIS  OF  THE  NECKLACE  ALGORITHM 
AND  ITS  APPLICATIONS 


Douglas  M.  Matty 


June  1999 


Thesis  Advisor: 
Second  Reader: 


Harold  M.  Fredricksen 
Craig  W.  Rasmussen 


Approved  for  public  release;  distribution  is  unlimited. 


X>nCQDAIJlTriirapBCTBD4 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  Including  the  time  for  reviewing  Instruction,  searching 
existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding 
thlo  burden  estimate  or  any  other  aspect  of  this  collection  of  Information,  including  suggestions  for  reducing  this  burden,  to  Washington  headquarters 
Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of 
Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188)  Washington  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  b/a/ikj 

2.  REPORT  DATE 

June  1999 

3.  REPORT  TYPE  AND  DATES  COVERED 

Master's  Thesis 

4.  TITLE  AND  SUBTITLE 

ANALYSIS  OF  THE  NECKLACE  ALGORITHM  AND  ITS  APPLICATIONS 

5.  FUNDING  NUMBERS 

6.  AUTHOR{S) 

Matty,  Douglas  M. 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

8.  PERFORMING 

ORGANIZATION  REPORT 
NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Department  of  Mathematics,  Naval  Postgraduate  School,  Monterey 
93940 

CA 

10.  SPONSORING  AGENCY 
REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the 
official  policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government. 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (maximum  200  words) 

A  commonly  studied  problem  xn  the  field  of  cryptography  is  the  Discrete  Logarithm 
Problem.  This  problem  is  also  referred  to  as  the  "distance"  problem.  Basically,  one 
would  like  to  know  where  a  particular  binary  n-tuple  is  in  a  list  combining  all  of  them, 
represented  as  powers  of  some  primitive  element,  or  equivalently  what  is  the  distance 
between  a  given  pair  of  n-tuples  in  a  siinilar  representation.  A  de  Bruijn  sequence  is  a 
well-known  periodic  binary  sequence  in  which  every  n-tuple  from  0  to  2"-l  appears.  Our 
goal  is  to  better  understand  the  "prefer-ones"  de  Bruijn  sequence.  Ultimately,  we  wish 
to  understand  where  each  of  the  binary  n-tuples  appears  in  that  sequence.  Using  the 
Necklace  Algorithm,  the  sequence  of  n-tuples  can  be  generated.  This  list  has  some 
special  properties  that  allow  us  to  perform  the  required  analysis  to  locate  the  n-tuples 
by  an  association  into  classes.  We  partition  the  binary  n-tuples  into  necklace  classes 
according  to  the  longest  substring  of  ones  appearing  on  the  n-tuple.  We  then  count  how 
many  n-tuples  appear  in  the  sequence  for  the  first  time  as  members  of  a  necklace  class 
containing  no  longer  strings  of  ones. 


14.  SUBJECT  TERMS 

Cryptography,  Combinatorics,  Discrete  Mathematics,  Analysis  of  Algorithms 


15.  NUMBER  OF 
PAGES 

53 


17.  SECURITY  CLASSIFICATION  OF  18.  SECURITY  CLASSIFICATION  OF 
REPORT  THISPAGE 

Unclassified  Unclassified 


19.  SECURITY  CLASSIFI- CATION 
OF  ABSTRACT 

Unclassified 


16.  PRICE  CODE 


20.  LIMITATION 
OF  ABSTRACT 

UL 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


ii 


Approved  for  public  release;  distribution  is  unlimited 

ANALYSIS  OF  THE  NECKLACE  ALGORITHM 
AND  ITS  APPLICATIONS 


Douglas  M.  Matty 
Captain,  United  States  Army 
B.S.,  United  States  Military  Academy,  1990 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  APPLIED  MATHEMATICS 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
June  1999 


Michael  A.  Morgan,/^hairman 
Department  of  Mamematics 


111 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


iv 


ABSTRACT 


A  commonly  studied  problem  in  the  field  of  cryptography 
is  the  Discrete  Logarithm  Problem.  This  problem  is  also 
referred  to  as  the  "distance"  problem.  Basically/  one  would 
like  to  know  where  a  particular  binary  n-tuple  is  in  a  list 
combining  all  of  them,  represented  as  powers  of  some 
primitive  element,  or  equivalently  what  is  the  distance 
between  a  given  pair  of  n-tuples  in  a  i similar 
representation.  A  de  Bruijn  sequence  is  a  well-known 
periodic  binary  sequence  in  which  every  n-tuple  from  0  to 
2"-l  appears.  Our  goal  is  to  better  understand  the  "prefer- 
ones"  de  Bruijn  sequence.  Ultimately,  we  wish  to  understand 
where  each  of  the  binary  n-tuples  appears  in  that  sequence. 
Using  the  Necklace  Algorithm,  the  sequence  of  n-tuples  can 
be  generated.  This  list  has  some  special  properties  that 
allow  us  to  perform  the  required  analysis  to  locate  the 
n-tuples  by  an  association  into  classes.  We  partition  the 
binary  n-tuples  into  necklace  classes  according  to  the 
longest  substring  of  ones  appearing  on  the  n-tuple.  We  then 
count  how  many  n-tuples  appear  in  the  sequence  for  the  first 
time  as  members  of  a  necklace  class  containing  no  longer 
strings  of  ones. 
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I .  INTRODUCTION 


A.  BINARY  N-TUPLES 

The  set  of  binary  n-tuples  is  the  set  of  all  strings  of 
O's  and  I's  of  length  n.  Usually  these  are  ordered  from 
00...0  to  11...1.  These  n-tuples  can  be  found  on  necklaces 
defined  by  cyclic  rotation.  These  necklaces  are  Lyndon 
words  of  length  d,  where  d\n.  Example  for  w  =  4: 

(0)  (GOOD  (0011)  (01)  (0001)  (1) 

These  are  the  Lyndon  words  (0)  and  (1)  of  length  one,  (01) 
of  length  two,  and  (0001) , (0011)  and  (0001)  of  length  four. 

B.  NECKLACE  ALGORITHM 

1.  Definition  of  Necklace  Algorithm 

The  Necklace  Algorithm  is  an  algorithm  for  generating 
the  necklaces  of  beads  of  length  n  in  two  colors,  i.e.,  the 
binary  necklaces.  Fundamental  to  the  Necklace  Algorithm  is 
the  operation,  defined  as: 

0  ( ai ,  3.2/ ... /  3n)  =  (bi , ba, ...,  b«)  / 

where  b,-  =  a,-  for  i  =  and  j  is  determined  as  the 

largest  subscript  such  that  ay>0  and  a^  =  0  for  all  k>j.  Then 
by  =  ay-1,  and  by+f  =  b,  for  7  =1, 2, ..., w-j. 
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Necklace  Algorithm  (for  necklaces  of  length  n) . 

Step  0.  The  initial  necklace,  Vq  ,  is  11...1  =  1”; 

Step  1.  0-step:  Form  ©(F^”)  to  find  F”^.,  ; 

Step  2.  J-check:  The  resulting  string  is  the  next 
necklace  if  and  only  if  j\n. 

Step  3.  If  we  have  not  found  the  last  necklace,  0”, 

by  steps  1  and  2  above,  return  to  step  1.  [1] 

Figure  1.  Necklace  Algorithm 

Note:  If  step  2  is  omitted,  we  produce  all  prenecklaces 
(defined  later)  .  If  step  2  is  modified  such  that  j  =  n,  we 
produce  a  list  of  all  Lyndon  words  of  length  n. 

The  following  output  of  the  Necklace  Algorithm  using 
the  three  options  for  the  J-check  is  provided  for  w  =  4: 


No  y-check 

j\n 

7  =  w 

1111 

1111 

1110 

1110 

1110 

1101 

1100 

1100 

1100 

1010 

1010 

1001 

1000 

1000 

1000 

0000 

0000 

Table  1.  Necklace  Algorithm  Output  («  =  4) 

Thus  the  Necklace  Algorithm  produces  a  listing  of  the 
n-tuples  appearing  from  the  largest  to  the  smallest  and  with 
d  of  the  n-tuples  appearing  on  each  Lyndon  word  of  length  d. 
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Additionally,  the  Lyndon  words  that  fom  the  necklaces  can 
be  juxtaposed  in  order  to  create  the  de  Bruijn  "prefer- 
one's"  sequence  as  follows:  1^1110^1100-10^1000^0.  [2] 

2.  Properties  of  the  Necklace  Algorithm 
The  Necklace  Algorithm  output  is  the  list  of  n-long 
necklaces,  where  the  largest  n- tuple  on  each  necklace 
represents  the  necklace  and  all  the  necklaces  appear  in 
lexicographic  order  from  largest  to  smallest.  Other  authors 
order  necklaces  and  their  respective  Lyndon  words  from 
smallest  to  largest  and  represent  a  necklace  by  its  least 
representative.  The  results  are  equivalent  and  it  is  easy 
to  translate  between  them.  Inherent  in  this  list,  the 
n-tuples  are  partitioned  into  classes  defined  according  to 
the  length  of  the  initial  string  of  I's  in  the  n- tuple.  We 
make  some  clarifying  definitions. 

Definition  1.  An  n-tuple  of  length  n  has  length 


parameter  /, 


where  /  =  >^ 


n 

2 

n-\ 


n  is  even 
n  is  odd 


and  is  a  member  of  the 


class,  m.  These  parameters  also  indicate  the  number  of 
leading  I's  as  there  are  /-w  leading  l^s  before  the  first  0. 
Each  class  then  has  the  general  form  with  the  respective 
parameters  as  follows: 
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DOF 


1...10a,_„+2  such  that  l-n  <  m  <1. 

Note:  The  extreme  special  cases  are  the  n-tuple,  1”, 
where  m  =  l-n  since  there  is  no  zero  following  the  lead 
string  of  ones,  and  0",  where  m  -  I  since  no  ones  preceding 
the  first  zero. 

l-m 

,.The  leading  substring  is  defined  as  and  is 

common  to  all  members  of  a  respective  m  class  for  a  given  /. 
To  consider  the  number  of  members  of  a  class,  we  also 
consider  the  degrees  of  freedom  {DOF),  i.e.,  those  bits 
that  are  unspecified  by  the  characteristics  of  the  class. 
We  consider  a  particular  n-tuple  class  as  an  example: 

V  —  1 1 1  OCX5(X60t70t80t9(XloOtll 
7:  5 

«  (Odd)  :  11  (e.g.  2/+1  =  2(5)+l  =  11) 

m:  2  (e.g.  I- [l-m)  =  5-3  =  2) 

DOF:  1  (e.g.  2/+1- [  (/-m) +1]  =  ll-[(5-2)+l]  =  7) 

A  subtle  point  is  to  be  noted.  The  leading  I's  in  a 
class  m  will  increase  in  number  with  7,  and  thus  with  n, 
while  m  remains  fixed.  So  an  11 -tuple  in  class  3  has  two 
leading  I's,  while  a  13-tuple  in  class  3  has  three  leading 


"I's."  The  ll-tuple  has  8  degrees  of  freedom  while  the 
13 -tuple  has  9  degrees  of  freedom. 

The  general  form  of  the  output  pf  the  Necklace 
Algorithm  is  refined  to  represent  the  output  into  the 
respective  m  classes.  The  first  output  in  class  m  appears 
by  applying  ©  to  the  last  output  of  the  class  m-1.  The 
last  element  of  the  class  m-1  is  denoted  as  [/-/m+1]  and  the 
first  element  of  the  class  m  is  then  denoted  as  0[/-w+l], 
where 


0[/-m+l]  =  ll... 1011... 1011. ..10...  111...  , 

and  the  last  output  of  the  Necklace  Algorithm  for  the  class, 
m,  referenced  as  [l-m] ,  has  the  form: 

l-m 

111000...0  . 

3.  Prenecklace  Generation 

A  prenecklace  appearing  by  step  1  of  the  Necklace 
Algorithm  is  defined  as  follows: 

Definition  2.  An  n-vector,  V" ,  is  a  prenecklace'  if 
there  exists  some  Ar-vector  F*  such  that  =  f(F",  F*) ,  where 
F”'"*  is  a  necklace  of  length  n+k.  (Note:  f()  denotes  the 

operation  of  concatenation  between  F”  and  F*)  [3]  . 
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Theorem  1 . 


The  vector  F”  =  ViV2...v„  is  a  prenecklace  if 
and  only  if  Vi...  v„_a+i  >  Vi.„v„  for  all  k  such  that  l<^<w.  [3] 

The  Lyndon  words  of  length  d  for  d\n,  when  extended  to 
length  n,  form  all  the  necklaces  of  length  w.  The  Lyndon 
words  of  length  d  for  d  <  n,  when  extended  to  length  n,  form 
all  the  prenecklaces  of  length  n.  [4]  We  describe  the 
relationship  between  the  prenecklaces,  necklaces  and  the 
Lyndon  words  in  the  Table  2  and  Figure  2  below  for  «  =  4. 


Prenecklaces 

Lyndon  Words 

Necklaces 

1111 

1 

1 

1110 

1110 

1110 

1101 

110 

1100 

1100 

1100 

1010 

10 

10 

1001 

100 

1000 

1000 

1000 

0000 

0 

0 

Table  2.  Comparison  of  Type  of  Words 


There  is  a  bijective  mapping  of  the  set  of  all  Lyndon 
words  of  lengths  one  through  n  onto  the  set  of  all 
prenecklaces  of  length  n  by  taking  each  of  the  Lyndon  words 
of  length,  d  <  n  and  extending  them  to  length  n.  These  are 
then  the  outputs,  as  defined  by  0,  to  make  the  «-long 
prenecklaces,  [4] 
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Ruskey  shows  that  to  count  P2(p)  /  the  binary 
prenecklaces  of  length  n,  we  need  only  sum  L^if)  where 
\<i<n,  all  of  the  Lyndon  Words  of  smaller  or  equal  length: 

P,(r,)  =  Y,hO) .  [31 

i=l 

Equation  1 

Likewise  there  is  a  bij  action  from  the  set  of  all 
Lyndon  words  of  length  d,  such  that  d\n,  onto  the  necklaces 
of  length  n.  To  count  N^in) ,  the  necklaces  of  length  n,  we 
need  only  to  sum  all  of  the  Lyndon  words  of  length  d  such 
that  d\n: 

[51 

Equation  2 

These  relationships  imply  that  the  necklaces  are  a 
subset  of  the  prenecklaces.  This  is  illustrated  in  the 
following  figure: 
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Figure  2.  Necklace  Algorithm's  Reduction  of  Search  Space 

Prior  to  the  publication  of  the  Necklace  Algorithm  [1], 
generation  of  the  necklaces  was  accomplished  by  selecting  an 
n-tuple  and  performing  n  cyclic  shifts  to  determine  the 
representative  of  the  necklace  class.  Graphically,  this 
would  be  similar  to  picking  a  point  in  the  space  F” ,  and 
moving  from  point  to  point  n  times  until  one  point  was 
determined  that  is  the  representative  of  the  necklace  class. 
If .  By  step  0,  the  Necklace  Algorithm  begins  the  search  for 
necklaces  in  iV”.  By  0,  the  search  for  the  next  necklace 
continues  in  the  set  of  prenecklaces,  P”.  [3]  By 
constraining  our  search  to  this  subset  of  all  n- tuples,  we 
realize  a  savings  in  the  number  of  strings  that  are 
considered  as  a  possible  necklace.  Then  the  J-check  removes 
the  need  to  perform  the  cyclic  rotation  to  test  whether  the 
n-tuple  represents  the  necklace. 
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Table  3  demonstrates  two  alternative  methods  of 
counting  the  prenecklaces  of  length  n.  First  the  partition 
based  on  the  niomber  of  Lyndon  words  of  length  i,  L2{i) ,  is 
given.  [6]  Meanwhile,  the  Necklace  Algorithm  partitions  the 
prenecklaces  into  their  respective  m  classes.  Table  3  is 
provided  for  the  case  k  =  11. 


i 

L2  U) 

m 

10(11, AW)  1 

- 

- 

-6 

1 

1  . 

2 

-5 

1 

2 

1 

-4 

2 

3 

2 

-3 

4 

4 

3 

-2 

8 

5 

6 

-1 

16 

6 

9 

0 

32 

7 

18 

1 

61 

8 

30 

2 

105 

9 

56 

3 

128 

10 

99 

4 

53 

11 

186 

5 

1 

^2(11) 

412 

10(11)1 

412 

Table  3.  Comparison  of  Counting  Prenecklaces 


A  particular  prenecklace  will  appear  in  different 

"\  . .  -i'- 

locations  in  the  two  tables,  (e.g.,  1  is  located  in  the 

rows  /  =  1  and  m  =  -6,  while  0^^  appears  in  row  i  -  1  and 
m  =  5  .)  The  discussion  of  how  to  count  the  prenecklaces 
with  respect  to  the  class  m  appears  later  and  is  the  major 
contribution  of  this  thesis.  By  carefully  counting  the 
savings  in  the  Necklace  Algorithm  listing  we  determine  the 
location  of  any  particular  n-tuple  in  the  prefer-ones  de 
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Bruijn  sequence.  The  savings  measure  for  each  m  class  and 
its  difference  from  the  power  of  2  described  by  the 

respective  degrees  of  freedom  is  the  number  of  necklaces 
appearing  as  enumerated  in  the  second  part  of  Table  3. 

4 .  Necklaces 

By  step  2,  the  J-check  removes  those  prenecklaces  that 
are  not  necklaces.  From  our  statements  in  the  previous 
section  (section  3),  step  2  therefore  removes  the  Lyndon 
words  of  length  i\n.  This  results  in  additional  savings  to 
those  observed  for  step  1.  The  following  theorem  verifies 
the  claim  that  step  2  produces  only  necklaces: 

Theorem  2:  A  vector,  V ,  satisfying  step  1  of  the 

Necklace  Algorithm,  is  a  necklace  if  and  only  if  j\n. 

Proof:  Given  the  vector,  F”  =  Vi  V2...v„,  suppose  n=  tj. 

then  vi  V2...v„  =  (Vi...Vy-iO) Clearly  this  is  a  necklace  only 
when  Vj . .  .Vy_,  0  >  v,. . .  OVj . . . v,._j  V/  9^  0  • 

Now  let  n  =  tj+s  with  1<5<7,  and  assume  without  loss  of 

generality  that  v,  >  v,.  ...v^._,0v,  as  before. 

Then  v,  ...v^v,  >  v,  ...Vj._,0v,  ...v,  /  since  dropping  the  first  5 

bits  in  each  string  yields  v,  ...v^._,0v,  ...v^  and  the 

prenecklace  is  not  a  necklace. □ 
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Alternative  to  Equation  2  above,  we  also  have 
MacMahon' s  formula  for  the  number  N(n, 2)  of  binary  necklaces; 

.  '  [7] 

^d\n 

Equation  3 

Here  is  Euler's  totient  function  and  the  summation  is 

over  all  divisors  of  the  necklace  length,  n. 

As  mentioned  in  t'he  discussion  of  prenecklaces  above, 
the  formula  does  not  partition  the  necklaces  into  classes  as 
does  the  Necklace  Algorithm.  However,  the  Necklace 
Algorithm  does  not  emamerate  the  necklaces  in  classes  as  we 
would  like.  (Note:  see  Appendix  A) 

We  perform  some  elementary  analysis  to  compute  the  size 
of  the  classes  for  this  parameter.  From  the  definition  of 
the  classes,  evidently  only  one  necklace  class  has  a  1  in 
the  last  position,  namely  1".  Any  other  prenecklace  that 
ends  in  a  1  can  simply  be  cyclically  rotated  to  the  right  to 
place  the  1  in  the  first  position,  making  the  resulting 
vector  bigger  than  the  original  vector,  e.g.,  1101  -»  1110 
and  1110  >1101. 

From  the  previous  statements,  an  elementary  upper  bound 
for  the  number  of  members  in  a  class  can  be  determined. 
Since  the  last  position  of  a  necklace  must  be  a  0,  an  upper- 
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bound  is  The  savings  described  above  shows  how  far 
below  this  power  of  2  the  number  of  necklaces  actually  is. 
By  careful  analysis  of  the  Necklace  Algorithm,  this  upper 
bound  can  be  tightened. 


II.  COUNTING  THE  NECKLACE  ALGORITHM 

A.  PRENECKLACE  COUNT 

1.  Table  of  Data  for  Missed  n-tuples 

Our  ultimate  goal  is  to  better  understand  the  “prefer- 
ones”  de  Bruijn  sequence.  We  wish  to  show  where  each  of  the 
binary  n-tuples  appears  in  the  sequence.  The  first  step  is 
to  determine  to  which  necklace  the  n- tuple  belongs.  For 
each  necklace  it  is  possible  to  determine  its  class,  as  the 
binary  n-tuples  have  been  partitioned  into  classes  according 
to  the  longest  leading  substring  of  ones  by  the  Necklace 
Algorithm.  If  we  can  count  how  many  necklaces  appear  in  the 
sequence  for  the  first  time  as  members  of  a  necklace  class 

containing  no  longer  strings  of  ones  we  will  have  found  the 

location  of  the  n- tuple.  We  choose,  in  fact,  to  count  how 

many  necklaces  fail  to  appear  in  this  way  (i.e.,  those 
necklaces  that  appear  earlier  in  the  sequence  as  members  of 
other  necklace  classes  containing  longer  strings  of  ones  as 
well).  That  is,  our  analysis  proceeds  as  we  count  the 

strings  that  are  missing  from  the  /w-class  of  size 

(Note:  see  Table  3  and  subsequent  discussion) . 


The 

computer  code  used  to 

implement 

the 

Necklace 

Algorithm 

is  modified  to  count 

the  number 

of 

n-tuples 
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omitted  (in  the  above  sense)  from  inclusion  in  the 
respective  m  class  in  the  list  of  prenecklaces.  For  the 
respective  m  classes  the  number  of  possible  n-tuples  with 
n- {l-m-1)  degrees  of  freedom  is  Not  all  of  these 

appear  employing  the  Necklace  Algorithm  as  we  have  noted. 
The  number  of  omissions  is  used  as  a  performance  measure  for 
the  efficiency  of  the  algorithm  as  stated  above  and  in 
Ruskey.  [3]  The  entire  table  of  data  is  included  in 
Appendix  A.  These  tables  in  Appendix  A,  aS  presented,  are 
segregated  by  odd  and  even  values  for  the  length  of  the 
n-tuples.^  Our  analysis  will  be  for  the  case  in  which  n  is 
odd.  When  n  is  even  the  analysis  is  similar.  The  data  in 
Table  4  and  Table  5  is  extracted  from  Appendix  A.  This  data 
specifies  the  number  of  n-tuples  missing  from  the  list  of 
prenecklaces  of  length  n  in  the  class  m. 


n  ^  2l  +  1 

l-m 

1 

9 

11 

13 

15 

/-4 

0 

0 

0 

0 

0 

/-I 

3 

3 

3 

3 

3 

1-2 

21 

23 

23 

23 

23 

1-3 

63 

104 

128 

136 

138 

/-4 

— 

255 

459 

704 

1-5 

— 

-- 

1023 

1930 

2871 

1-6 

— 

— 

— 

4095 

7926 

1-1 

— 

— 

--  , 

— 

16383 

Table  4 .  Number  of  Missing  n-tuples  (Odd) 
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n  =  2  1 

l-m 

6 

8 

10 

12 

14 

l-O 

0 

0 

0 

0 

0 

1-1 

1 

1 

1 

1 

1 

1-2 

9 

9 

9 

9 

9 

1-3 

31 

48 

56 

58 

58 

l-A 

- . 

127 

221 

288 

314 

1-5 

— 

— 

511 

946 

1353 

1-6 

— 

— 

--  ■ 

2047 

3920 

1-1 

1  - 

— 

— 

— 

8191 

Table  5.  Number  of  Missing  n-tuples  (Even) 

Table  4  and  Table  5  provide  several  insights  for  our 
further  analysis.  First,  there  are  no  missing  n-tuples  when 
the  number  of  leading  ones  is  greater  than  or  equal  to  /, 
that  is,  when  /»  <  1.  Second,  when  there  are  no  leading 
ones,  i.e.,  m  =  I,  the  only  appearing  n-tuple  is  00. ..0  =  0". 
Thus,  there  are  n-tuples  missing  that  have  the  initial 

string  of  one's  being  zero  long.  For  example  for  n  =  6, 
/=3,  the  last  n-tuple  from  the  Necklace  Algorithm  belonging 
to  the  class  m  =  3  is  000000,  which  has  /-w  =  /-3  =  3-3  =  0 , 
i.e.,  the  leading  number  of  ones  is  zero  long,  the  degrees 
of  freedom  {DOF)  is  five  and  31  n-tuples  are  omitted  in 
this  class.  The  dashed  lines  in  the  tables  indicate  that  no 
value  is  defined  for  the  given  parameters  as  the  numbers  are 
not  meaningful.  A  third  observation  is  that  the  numbers  in 
the  table  seem  to  increase  to  a  value  and  remain  at  that 
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value.  We  call  this  the  steady  state  number  for  that  m 

class.  The  values  when  steady  state  has  been  reached  are 
indicated  in  bold.  This  steady  state  number  reflects  the 
maximum  savings  afforded  by  performing  the  Necklace 
Algorithm  for  that  class.  Implicitly,  the  first  sub-problem 
we  face  is  to  count  the  size  of  the  steady  state  set  of 

missing  n- tuples  for  a  given  class  m.  Further,  we  show  for 

which  value  of  n  the  steady  state  occurs  for  a  given  class 

m. 

The  notion  of  a  class  reaching  steady  state  is 
demonstrated  in  Table  6.  Consider  the  following  n-tuples 
from  the  class  /w  =  2,  which  are  missing  from  the  list  of 
prenecklaces  for  the  lengths  of  nx  7,  9,  and  11.  We  do  not 
include  the  value  of  n  =  5,  which  for  the  class  m  =  2  has 
zero  leading  I's.  Thus  there  are  15  missing  members  in  the 
class,  namely  the  non-zero  5-tuples  with  a  leading  0. 
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It: 

«  =  9 

K  =  11 

Count 

1011111 

110111111 

11101111111 

1 

1011110 

110111110 

11101111110 

2  . 

1011101 

110111101 

11101111101 

3 

1011100 

=> 

110111100 

11101111100 

4 

1011011 

=> 

110111011 

11101111011 

5 

•  => 

110110111 

11101101111 

.  6 

1011010 

=> 

110111010 

11101111010 

7 

1011001 

=> 

110111001 

11101111001 

8 

1011000 

'  => 

110111000 

11101111000 

9 

1010111 

=> 

110101111 

11101011111 

10 

1010110 

=> 

110101110 

11101011110 

11 

1010011 

=> 

110100111 

11101001111 

12 

1001111 

■=> 

110011111 

11100111111 

13  / 

1001110 

=> 

110011110 

11000111110 

14 

1001101 

llOOlllOl 

11100111101 

15  : 

=> 

110011011 

11100111011 

'16 

1001100 

110011100 

11100111100 

17  : 

1001011 

=> 

110010111 

11100101111 

18 

1001010 

110011010 

11100111010 

19 

1000111 

=> 

110001111 

11100011111 

20  : 

1000110 

=> 

110001110 

11100011110 

21 

1000101 

=> 

110001101 

11100011101 

22 

lOOOOll 

=> 

110000111 

11100001111 

23 

21 

23 

23 

Table  6. 

Increasing  Missing 

n-tuples  to  Steady  State 

At  n  = 

7  the 

m  -  2  class 

is  not  yet  at 

steady  state; 

>  is  the 

first 

value  of  n 

that  the  class 

»2  =  2  is  at 

steady  state.  At  n  =  11,  the  Necklace  Algorithm  produces  a 
list  of  necklaces  that  remains  at  steady  state  and  there  are 
twenty-three  n-tuples  omitted  for  the  class  m  =  2  for  each 
succeeding  value  of  n.  Comparing  the  three  columns,  we  see 
that  each  missing  n-tuple  can  be  mapped  {=>)  to  an  n- tuple 
in  the  successive  column  by  adding  a  1  to  the  leading 
substring  (this  is  required  for  to  maintain  the  class 
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structure) ,  and  adding  an  additional  1  to  the  longest  string 
of  I's  in  the  trailing  substring.  Proceeding  from  w  =  7  to 
n  =  9,  there  occasionally  are  two  opportunities  to  add  the 
second  1.  It  is  interesting  to  note  that  from  k  =  9  to 
n  =  \\  each  missing  n-tuple  maps  to  one  and  only  one  n-tuple 
in  the  next  column.  For  example,  1011011  produces  both 
110111011  and  11011011.  In  section  3  we  show  when  such 
possibilities  for  two  successors  end  and  steady  state  is 
therefore  achieved. 

2.  Plot  of  Missing  n-tuples  Relative  to  m 

An  approximation  to  the  maximum  number  of  missing 
n-tuples  is  made  by  analyzing  the  data  in  Appendix  A.  We 
extract  from  the  data  the  first  length  n  for  each  class  m 
at  which  steady  state  is  reached.  Our  analysis  includes  the 
values  for  the  parameter  n  (2/+1  for  odd  and  2/  for  even),  /, 
m,  l-m,  the  degrees  of  freedom  and  the  maximum  number  of 
missed  n-tuples.  The  number  of  n-tuples  missed  is 
transformed  using  the  logarithm  base  two  to  aid  in  the 
scaling  on  the  graph.  (Base  2  seems  a  natural  choice  since 
we  are  concerned  with  a  binary  alphabet.) 
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21+1 

1 

m 

l-m 

DOF 

#missed 

loga  0 

3 

1 

1 

0 

2 

3 

1.585 

9 

4 

2 

2 

6 

23 

4.524 

15 

7 

3 

4 

10 

138 

7.109 

K! 

10 

4 

6 

■D 

740 

9.531 

27 

13 

5 

8 

18 

3720 

11.861 

33 

16 

6 

10 

22 

17936 

14.131 

Table  7.  Steady  state  for  odd  values  of  n 


21 

1 

m 

l-m 

DOF 

ttmissed 

logaO 

6 

3 

2 

1 

4 

9 

3.170 

12 

6 

3 

3 

8 

58 

5.858 

18 

9 

4 

5 

12 

324 

8.340 

24 

12 

5 

7 

■EB 

1672 

10.707 

30 

15 

7 

8 

21 

8208 

13.003 

36 

18 

8 

10 

25 

38944 

15.249 

Table  8.  Steady  state  for  even  values  of  n 


There  is  evidently  a  strong  relationship  between  I  and 
the  nimiber  of  missing  n- tuples.  The  following  graphs  Of  / 
vs.  loga  (#missing  n- tuples)  illustrates  this  nearly  linear 
relationship.  For  reference,  a  line  of  the  form  y  =  x  is 
plotted  in  dashed  lines  as  well. 
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Log2(#missed) 


Odd  Steady  States 


14  7  10  13  16 


/ 


Figure  3 .  Regression  of  Steady  States  (Odd) 


Even  Steady  State 


3  6  9  12  15  18 


/ 


Figure  4 .  Regression  of  Steady  States  (Even) 
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This  linearity  is  checked  with  a  linear  regression  of 
the  data  where  the  expected  sum  of  squares  yields  a 
correlation  of  determination  (R^)  of  0.99790  for  the  odd 
length  n- tuples  and  0.9989  for  the  even  n-tuples.  This 

leads  us  to  propose  2^  «  steady  state  value.  However,  it  is 

also  evident  that  this  is  only  a  loose  upper  bound  as  /  (and 
therefore  n)  grows.  Our  further  analysis  will  also  serve  to 
tighten  this  bound. 

3.  Finding  Steady  State 

The  Necklace  Algorithm  proceeded  using  m  as  a 
parameter.  The  above  argument  establishes  /  as  a 
determining  factor  for  the  number  of  missing  n-tuples  (which 
is  also  in  a  linear  relationship  with  n) .  For  the  odd  values 
of  n,  steady  state  is  reached  for  .«  =  3,9,15,21,27,33,...  . 
Evidently  steady  state  is  reached  when  «s3  (modulo  6).  '  :  We 
consider  the  structure  of  the  n-tuples  when  the  steady  state 
is  reached.  The  predecessor  and  the  first  n-tuple  in  the 
class  where  steady  state  is  reached  exhibit  a  similar 
property  among  the  classes. 


21 


m 

/ 

2/+1 

[/-/w+1] 

■ 

eu-m+i] 

1 

fl 

3 

i 

0(100) 

1-1 

'’ooo 

2 

4 

9 

M 

©(moooooo) 

1-2 

noiioiio 

3 

7 

15 

1-2 

©(iTmoooooooooo) 

/-3 

imoiiiioiiiio 

4 

10 

21 

1-3 

0(i  1 1 1 1 1  ioooooooooooooo) 

1-4 

mirioiii  11 101111110 

Table  9.  Examples  of  First  Member  of  Steady  State  Class 


By  the  Necklace  Algorithm,  the  last  one  in  an  n-tuple 
is  complemented  producing  /-m-1  ones  and  a  zero.  It  is 
apparent  that  the ■ length,  l-m,  of  this  substring  pattern 
divides  n  with  a  cofactor  of  3.  We  also  consider  larger 
values  for  n  to  ensure  that  steady  state  is  maintained. 
Table  10  illustrates  the  first  output  of  the  Necklace 
Algorithm  for  the  class  m  -  2  for  various  lengths  of  n. 


n 

/ 

[/-/w+1] 

■ 

0[/-/w+l] 

7 

3 

/-I 

0(nooooo) 

1-2 

1010101 

9 

4 

1-1 

©(moooooo) 

1-2 

110110110 

11 

5 

i-i 

©(imooooooo) 

/-2 

11101110111 

13 

6 

/-I 

©(nmoooooooo) 

1-2 

imoiiiiom 

Table  10.  Exan^les  of  First  Member  of  Class  m  =  2 
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It  can  be  confirmed  from  the  data  in  Table  A,  for  the 
class  1-2,  steady  state  is  reached  at  2/+1  =  9/  and  continues 


for  all  larger  values  of  /.  The  first  output  of 


the 


Necklace  Algorithm  for  the  class  m  =  2,  and  for  the  value 
of  n  for  which  steady  state  is  reached  has  the  following 
form  {R  is  used  to  abbreviate  the  remainder) ; 


t-2  1-2  R 

•  rTlOMOna^  . 

The  term  remainder  indicates  that  part  of  the  leading 
substring  that  is  not  completely  copied.  Evidently,  this 
shows  that  if  /  is  large  enough  to  allow  the  leading 
substring  to  be  copied  two  '  times  and  leave  a  remainder  of 
length  three  or  greater  then  steady  state  will  be  reached 
for  the  class  m  =  2.  The  classes  of  larger  values  of  m 
also  exhibit  a  distinct  relationship  with  /  and  m.  For 

n  =  2/+l  =  31,  and  the  classes  of  m  =  1,2,..., 6  the  following 
similar  characteristics  appear  (see  Table  11) . 
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m 

/ 

Last  output 
class  /-I 

■ 

First  output 
class  1-2 

Length  of 
Remainder 

1 

15 

1 

0(1...  100000) 

1 

2 

15 

7-1 

0(000000) 

3 

3 

15 

5 

4 

15 

7-3 

0(l?7iooooo) 

7-4  7-4 

Oorrioiiiiiii 

7 

5 

! 

15 

i-4 

0(OlOOOOO) 

7-5  7-5 

OolTioiiiiiii 

9 

6 

15 

7-5 

0(Ciooooo) 

1 

Table  11.  Examples  of  First  Output  for  Classes  (/  =  15) 


(Note:  for  m  =  6,  steady  state  has  not  been  reached  for 
71  =  31 . ) 

This  suggests  the  following  theorem; 

Theorem  3 : 

Steady  state  for  the  0-step  of  The  Necklace  Algorithm 
will  be  reached  for  the  class  of  tw  at  77  =  2/+1  ,  when 

I  +  2> 3m  . 

Proof :  From  Table  10,  we  see  that  the  substring  "110" 
is  copied  three  times  for  2/  +  1  =  9.  For  larger  n  and  the 
same  m,  the  initial  string  l^'^O  will  only  be  copied  twice 
plus  the  first  2777-1  bits.  This  establishes  the  following 
inequality: 
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3(/-7M  +  1)  >  2/  +  1  =  « 

3/  —  'int  +  3>2/  +  l  =  w 
/  +  2  >  3m 

To  illustrate  Theorem  3,  we  verify  the  result  for  the 

n-tuple  of  length  31,  showing  that  the  class  of  m  -  6  and 

■  ■  ‘ 

the  number  of  leading  ones  (/-m  =  15-6)  is  not  at  steady 
state : 


Table  12.  Verification  of  Theorem  3 

(Note:  Steady  state  is  achieved  at  2/+1  =  33.) 

By  the  theorem,  steady  state  is  obtainable.  Because 
this  is  not  a  priori  apparent  we  note  this  specifically  in 
the  following  corollaries. 

Corollary  1 : 

For  a  given  class,  m,  the  Necklace  Algorithm  will  reach 
steady  state  for  some  I  >  L. 

Proof:  m  is  a  fixed  quantity.  Suppose  that  /  +  2<3m. 

Let  Z  =  3m.  Increase  /  such  that  /  >  Z.  □. 
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Corollary  2: 

Steady  State  is  reached  for  class  m  if  and  only  if: 

==  2  with  a  remainder  of  2m-l. 

It  is  helpful  to  consider  an  example  of  Corollary  2. 
Let  /  =  15,  then  n  =  31.  From  Table  12,  we  can  compare  the 
results  for  the  classes  m  =  1,2, 3, 4, 5, 6  using  the 

corollary. 


m 

2m-l 

\y(l-m  +  \)_ 

(/-w+1 )  mod  («) 

1 

1 

2 

1 

2 

.  •  3 

2 

3 

3 

5 

2 

5 

4 

7 

2 

■  7 

5 

9 

2 

9 

6 

11 

3 

1 

Table  13.  Coir^sare  Results  of  m  and  Steady  State 


(Note:  Steady  state  is  reached  for  ill  cases  but  m  =  6.) 

B.  NECKLACE  COUNT 

1 .  Table  of  Data  for  Omitted  Necklaces  and  n-tuples 

For  the  "savings"  realized  by  the  J-check,  we  first 
examine  the  difference  between  the  prenecklaces  and  the 
necklaces  of  length  n.  Equation  1  and  Equation  2 
respectively  compute  the  values  of  P"  and  if,'  respectively. 
These  values  are  presented  in  Table  14. 


% 


7-/W  +  1) 
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n 

pn 

N" 

P"-lf 

1 

2 

2 

0 

2 

3 

3 

0 

3 

5 

4 

1 

4 

8 

6 

2 

5 

14 

8 

6 

6 

23 

14 

9 

7 

41 

20 

21 

8 

71 

36 

35 

9 

127 

60 

67 

Table  14.  Prenecklaces ^  Necklaces ,  and  the  Difference 


The  number  P”  is  approximately  twice  the  number  If.  .But 
more  information  can  be  obtained  by  comparing  the  number  of 
prenecklaces  and  necklaces  for  each  of  the  w-classes  for  a 
fixed  Tt.  In  particular,  steady  state  becomes  apparent 
again.  Table  •  15  gives  values  partitioned  by  w-class  for 
the  case  «  =  31 . 
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m 

P(ni) 

N(m) 

P(m)-N(m) 

-16 

1 

1 

0 

-15 

1 

1 

0 

-14 

2 

1 

1 

-13 

4 

2 

2 

-12 

8 

4 

4 

-11 

16 

8 

8 

-10 

32 

16 

16 

-9 

64 

32 

32 

-8 

128 

64 

64 

-7 

256 

128 

128 

-6 

512 

256 

256 

-5 

1024 

512 

512 

-4 

2048 

1024 

1024 

-3 

4096 

2048 

2048 

-2 

8192 

4096 

4096 

-1 

16384 

8192 

8192 

0 

32768 

16384 

16384 

1 

65533 

32766 

32767 

2 

131049 

65522 

65527 

3 

262006 

130992 

131014 

4 

523548 

261728 

261820 

5 

1044856 

522240 

522616 

6 

2079218 

1038850 

1040368 

7 

4110418 

2052174 

2058244 

8 

8005876 

3990960 

4014916 

9 

15076085 

7492056 

7584029 

10 

26279377 

12974378 

13304999 

11 

38218348 

18604412 

19613936 

12 

35857846 

16913378 

18944468 

13 

11532736 

5064334 

6468402 

14 

269684 

97108 

172576 

15 

1 

1 

0 

Table  15.  Count  of  Step  Outputs  («  =  31) 


It  apparent  that  P(m)-N(m)  grows  by  a  factor  of  2  until 
m  >  1.  Then  when  the  Necklace  Algorithm  is  performed,  a 
necklace  representative  is  found  without  performing  the  n 


rotations  and  comparisons.  But  the  cost  of  the  algorithm  is 
still  about  twice  the  number  of  necklaces  found,  since  the 
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prenecklaces  have  to  be  determined  first.  Reviewing  Figure 
2,  the  number  of  prenecklaces  that  are  not  necklaces, 
P(n)-N(n),  is  evidently  on  the  order  of  the  number  of 
necklaces  [1] .  Thus  the  work  to  determine  the  necklaces  is 
cut  down  from  the  k2"  rotations  and  comparisons  to 


0(Ar(«))  =  0(2%) 


as  a  function  of  utilizing  the  Lyndon  words. 


C .  DEFINING  STEADY  STATE 

Steady  state  has  been  shown  to  occur  for  particular 
m-classes  of  prenecklaces  in  Theorem  3.  In  the  previous 
section,  the  prenecklaces  are  determined  to  be  «-long 
extensions  of  the  Lyndon  words  of  length  d,  for  1  <  d  <  n. 
Similarly,  the  necklaces  are  determined  by  extending  to 
length  n  the  Lyndon  words  of  length  d,  where  d\n.  We  now 
consider  partitioning  the  Lyndon  words  into  w-classes  to 
verify  this  notion  of  steady  state. 

From  the  partitioning  into  7w-classes  of  the 
prenecklaces  for  a  given  n,  \^e  can  count  the  number  of 
prenecklaces  for  a  given  length  of  n  using  Equation  2. 
Since  the  prenecklaces  are  f/-long  Lyndon  words  extended  to 
length  n,  we  can  count  the  number  of  prenecklaces  in  a  given 
class  by  counting  the  number  of  r/-long  Lyndon  words  that 
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begin  with  the  substring  Thus  we  can  enumerate  the 

Lyndon  words  for  each  /w-class  and  for  any  value  of  n. 
Table  16  lists  the  number  of  prenecklaces  of  length  nine 
determined  by  Lyndon  words  of  length  d  <  9  and  partitioned 
by  the  number  l-m  of  leading  I's: 


l-m 

d 

9 

8 

7 

6 

5 

4 

3 

2 

1 

0 

1 

I — 1 

0 

0 

0 

0 

0 

0 

0 

0 

mm 

2 

0 

0 

0 

0 

0 

0 

0 

0 

-.1 

0 

3 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

4 

0 

0 

0 

0 

0 

0 

1 

\ — 1 

1 

0 

5 

0 

0 

0 

0 

0 

1 

1 

2 

0 

6 

0 

0 

0 

0 

1 

1 

2 

3 

2 

0 

7 

0 

0 

0 

1 

1 

2 

4 

6 

4 

0 

8 

0 

0 

1 

1 

2 

4 

7 

10 

5 

0 

9 

0 

1 

1 

2 

4 

8 

14 

18 

8 

0 

Table  16.  Partitioning  Lyndon  Words 


(Note:  The  decreasing  values  for  l-m  represent  the 
lexicographic  order  imposed  by  the  Necklace  Algorithm.  The 
Necklace  Algorithm  proceeds  by  columns,  from  left  to  right 
in  the  table,  producing  each  »-long  extended  Lyndon  word 
(prenecklace) . 

From  the  table  it  is  evident  that  the  number  of 
prenecklaces  in  the  class  /w  =  2  is  (1+1+2+3+6+10+18  =)  41. 

In  Table  4,  steady  state  for  the  class  w  =  2  is  23  missing 
n-tuples.  We  verify  that  2^^-P(n)  =  64  -  41  =  23.  The 
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reader  can  easily  check  that  steady  state  is  reached  for  the 
class  m  =  1  (i.e.,  missing  3  n- tuples  in  the  column 
l-m  =  3  )  . 

Not  only  can  it  be  determined  how  many  prenecklaces  are 
missing  for  a  particular  class,  but  it  can  be  determined  how 
many  are  missing  for  a  particular  value  of  d.  ;  By 
precalculating  the  size  of  the  set  of  missing  n-tuples  at 
steady  state  for  each  7w-class,  it  is  possible  to  precisely 
determine  the  position  for  a  given  necklace.  A  necklace  in 
a  class  for  which  steady  state  has  been  achieved  can  easily 
be  calculated  using  the  known  value  for  a  precise  position. 
For  other  classes,  not  at  steady  state,  this  information  can 
be  helpful  in  establishing  a  tight  upperbound  for  the  number 
of  necklaces  in  a  given  class  before  initiating  a  sequential 
search  in  .the  list  of  necklaces. 
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III.  CONCLUSIONS  AND  SUGGESTIONS  FOR  ADDITIONAL  RESEARCH 

The  Necklace  Algorithm  has  been  shown  to  be 
significantly  more  efficient  than  a  naive  algorithm 
producing  all  of  the  necklaces  of  a  given  length  in 
lexicographic  order.  The  algorithm  also  has  been  shown  to 
produce  all  of  the  Lyndon  words  of  lengths  one  through  n  in 
lexicographic  order  as  prenecklaces.  By  developing  the 
notion  of  steady  state,  it  is  shown  to  be  possible  to  find 
the  exact  position  of  the  necklace  representative  of  any 
particular  n- tuple  given  that  the  leading  substring  has  at 
least  /  ones.  If  the  representative  necklace  has  fewer 
leading  I's,  a  table  is  available  to  count  the  missing 
necklaces  for  the  given  class  m.  The  class  must  be  at 
steady  state  to  ensure  that  the  table  size  is  feasible.  iThe 
condition  for  steady  state  is  shown  to  be  l  +  2<2m\ 
Additional  methods  must  be  employed  when  steady  state  has 
not  been  achieved.  ' 

To  demonstrate  the  application  of  our  results,  an 
example  is  given.  We  choose  an  n-tuple: 

=  OOOllOlllllOlOl. 

First,  we  determine  the  n-tuple^ s  parameters  using 
Definition'!.  The  longest  substring  of  I's  is  11111;  and 


33 


since  I  -  1  find  m  =  2.  Now,  we  check  for  the  steady 
state  condition: 

7+2  >  3(2)  — >  9  >  6  Class  m  =  2  is  at  steady  state. 

From  Appendix  A,  the  steady  state  values  for  odd  values 
of  n  are  3,  23,  138  in  the  classes  m  =  1,  2,  and  3.  The 
representative  of  the  necklace  class  is:  111110101000110. 

This  representative  has  no  internal  periodicity  so  the 

! 

J- check  would  have  shown  that  this  n- tuple  is  a  Lyndon  word 
of  length  =  15.  There  are  242  Lyndon  words  of  length  15 
and  with  l-m  =  5.  To  locate  the  exact  necklace  then  we 
would  enter  the  string  [/-m-1]  (  =  111111000000000).  We 
then  apply  the  Necklace  Algorithm  and  count  how  many 
necklaces  are  produced,  until  the  correct  one  is  found. 

From  this  example  and  our  discussion  in  the  previous 
section,  it  is  clear  that  a  better  understanding  of  the 
Lyndon  words  and  the  development  of  an  efficient  counting 
technique  for  a  partitioning  by  class  would  result  in  an 
improved  solution  for  this  version  of  the  distance  problem. 
With  the  knowledge  of  the  location  of  a  necklace  relative  to 
an  initial  position  in  the  "prefer-ones"  de  Bruijn  sequence 
a  distance  can  be  determined  [5]  and  refined  quickly  in  such 
a  manner  that  it  would  be  computationally  efficient.  [8] 
This  could  result  in  the  development  of  a  new  cryptographic 
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scheme.  Our  methods  and  results  may  also  have  application 
in  data  compression.  Certainly/  it  can  be  said  the  solution 
to  this  problem  is  valuable  for  both  academic  and  practical 
reasons. 
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APPENDIX  A.  TABLE  OF  MISSING  N-TUPLES  COUNTS 


mm 

mm 

'Msm 

mm 

15 

17 

19 

21 

23 

25 

27 

29 

31 

/-o 

HI 

■3 

■El 

■D 

0 

[mB] 

mHH] 

0 

imH] 

0 

WE 

wm 

iHKi 

3 

3 

3 

3 

3 

3 

3 

3 

3 

n 

ms 

WE 

WE 

23 

23 

23 

23 

23 

23 

23 

23 

23 

a 

Wigl 

wniH 

138 

138 

138 

138 

138 

138 

138 

138 

sm 

255 

459 

628 

704 

730 

738 

740 

740 

ism 

nigg 

2871 

3392 

3606 

3684 

3710 

3718 

3720 

asm 

g!iEH 

15631 

17038 

17606 

17822 

17900 

17926 

urn 

53717 

69872 

78079 

81646 

83100 

83670 

83886 

65535 

225253 

305619 

349970 

370379 

379118 

382732 

mm 

■n 

HI 

mil 

imi 

H^m 

^mH 

262143 

932995 

1315402 

1542563 

1652306 

1701131 

BIM 

o 

Ja. 

00 

2089242 

3830876 

5592057 

6710766 

7275055 

HW 

im 

mH 

mu 

imHi 

mHi 

imHH 

8369683 

15630420 

23543815 

28890516 

16777215 

33508823 

63477929 

WkW 

67108863 

134107119 

256902720 

BM 

268435455 

536601228 

BIM 

1073741823 

Table  17.  Missing  n-tuples  (nodd) 


7 

9 

KEB 

KEI 

15 

17 

mm 

21 

23 

25 

27 

29 

31 

GQII 

mu 

HKI 

mm 

HHEl 

imiHE] 

0 

0 

0 

■3 

mm 

mu 

IHHE] 

0 

HHE] 

IHHEI 

0 

mmKi 

^mHH 

0 

0 

mm 

m 

wm 

■Q 

0 

0 

_ ^ 

HHKl 

Hll^E] 

^^mm 

0 

0 

mm 

■Q 

H 

8 

2 

0 

H^^El 

0 

mmmj] 

mHHKi 

ism 

EQ 

76 

26 

2 

0 

0 

hhhq 

907 

941 

521 

214 

78 

26 

8 

2 

0 

0 

mm 

1407 

568 

216 

78 

26 

8 

mm 

15774 

21560 

16155 

8207 

3567 

1454 

216 

64118 

95600 

80366 

44351 

20409 

3614 

mi 

mH 

im 

i^^m 

mHi 

412040 

382407 

227161 

109743 

MO 

1040667 

1741634 

1761181 

1118709 

WIBI 

4175380 

7260737 

5346701 

BgM 

16731608 

29969106 

34881953 

BM 

66998256 

122795601 

Fgt-1 

I 

268165773 

IWM 

mi 

IH 

^m 

i^m 

HHl 

HHH 

HHl 

HHH 

I^^Hl 

HHHl 

HH^^I 

Table  18.  Forward  Difference  (« odd) 
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r~ 

L6 

IB 

IB 

B 

14 

16 

18 

20 

22 

24 

26 

_ 28 _ 

30 

\m 

WE 

BHU 

^m] 

0 

0 

0 

0 

BOB 

1 

wm 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

mm 

m 

K] 

mm 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

SB! 

El 

B 

mm 

58 

58 

58 

58 

58 

58 

58 

58 

58 

am 

lEH 

wfa 

288 

314 

322 

324 

324 

324 

324 

324 

324 

324 

ED 

ms 

1353 

1558 

1636 

1672 

1672 

1672 

mki 

6038 

7319 

00 

1— • 

VO 

00 

8208 

fsm 

8191 

15983 

26075 

33160 

36599 

38044 

38614 

38830 

38908 

32767 

110209 

146483 

177932 

179388 

ESB 

131071 

259975 

459035 

635182 

736127 

783764 

804579 

RCT 

524287 

1043452 

1892322 

2715891 

3222158 

3472217 

UTM 

7743350 

11486697 

mmsernm 

VSfM 

mu 

Him 

mHi 

mHi 

imn 

HHHf 

16747871 

31514779 

48164396 

33554431 

67037905 

127749326 

QQi 

134217727 

QSi 

5368709111 

Table  19.  Missing  n-tuples  (« even) 


— 

6 

rr- 

IB 

IB 

14 

16 

18 

22 

24 

26 

28 

30 

IH] 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

QB 

Bi] 

0 

[HE] 

miE] 

0 

imE] 

0 

mHK] 

0 

IH^HHi] 

0 

B 

H] 

0 

HE] 

0 

H^HI] 

0 

mimi] 

IHIH] 

0 

0 

HHHK!] 

0 

B 

■s 

8 

2 

i[miQ 

EHH] 

0 

0 

0 

0 

B 

67 

26 

8 

2 

0 

0 

0 

0 

0 

0 

B 

435 

407 

205 

78 

26 

8 

2 

0 

0 

B 

1873 

2118 

1281 

559 

216 

78 

8 

2 

B 

7792 

570 

216 

78 

B 

31842 

45600 

36274 

19233 

8611 

3605 

1456 

B 

128904 

199060 

176147 

100945 

47637 

20815 

I/-10 

519165 

823569 

250059 

Bf 

H 

HHf 

^m 

^mH 

HIH 

^HH 

2453799 

BBB 

HHI 

i^m 

HHH 

HHI 

HBHH 

HHHi 

8359264 

14766908 

16649617 

BEDi 

33483474 

60711421 

M4 

134045153 

QBI 

Table  20.  Forward  Difference  (« even) 


APPENDIX  B.  CODE  FOR  NECKLACE  ALGORITHM 

The  following  code  is  a  Java  program  that  implements 

the  Necklace  Algorithm: 

import  java. util.*; 
import  java.lang.*; 
import  java.io.*; 
public  class  Theta  { 

public  static  void  main  (String  args[])  { 
int  length  =31; 
int  k  =  1; 

int  sequence  []  =  new  int [length]; 

int  j  =  length-1; 

int  counter  =0; 

int  numNeck  =0; 

int  numType  =0; 

int  flag  =1; 

int  group  []  =  new  int [length+1] ; 
int  sum  =0; 

System. out. println ("Finding  all  binary  necklaces  of  length  ”+ 
length 

//  Sets  the  necklace  of  all  l*s. 
for  (int  a  =  0;  a  <  length;  a++)  { 
sequence [a] =k; 

System,  out .print (sequence [a] ) ; 

) 

System. out .println ("  ") ; 

//  Find  last  non  zero  element  and  decrement  it  : 

while  (sequence [0]  !=  0)  { 

j  =  length-1; 
while  (sequence[j]  ==  0) { 

j=j-l; 

■  '  . 

//  Change  last  one-bit  and  fill  register. 
sequence[j]  =  sequence [j ] -1; 
if  (j  !=  length-1)  { 

for  (int  i  =  counter+1;  i  <  length- j ;  i++) { 
sequence [j+i]  =  sequence [i-1] ; 

} 

if  (length  %{j+l)  ==  0)  {  //  J-check 

for  (int  a  =  counter;  a  <  length;  a++)  { 

,  System,  out .print (sequence [a] ) ; 

} 

■  } 

} 

}  ■■  ■  ■ 

} 
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APPENDIX  C.  TABLE  OF  LYNDON  WORD  COUNTED  BY  PARTITIONS 


d 

14 

13 

9 

8 

7 

6 

5 

4 

3 

2 

Bi 

m 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

m 

■5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

2 

i 

6 

0 

0 

0 

0 

0 

0 

0 

1 

1 

2 

3 

2 

0 

7 

0 

0 

0 

0 

0 

0 

1 

1 

2 

BBS 

4 

0 

$ 

0 

0 

0 

0 

.  0 

0 

0 

1 

1 

2 

4 

7 

10 

5 

0 

9 

0 

0 

0 

1 

1 

2 

4 

8 

14 

8 

0 

la 

0 

0 

0 

0 

1 

1 

2 

4 

8 

15 

26 

0 

11 

0 

0 

1 

1 

2 

4 

8 

16 

30 

50 

56 

18 

0 

12 

0 

0 

0 

1 

1 

2 

4 

.8 

16 

31 

58 

93 

96 

25 

0 

B 

0 

0 

1 

1 

2 

4 

8 

16 

32 

62 

0 

B 

1 

1 

2 

8 

16 

32 

63 

122 

221 

334 

299 

58 

0 

B 

1 

1 

2 

4 

8 

16 

32 

64 

0 

16 

1 

2 

4 

8 

16 

32 

64 

127 

477 

1194 

929 

135 

0 

17 

2 

4 

8 

32 

64 

128 

254 

B 

4 

8 

64 

128 

255 

506 

2893 

B 

B 

8 

16 

32 

64 

128 

256 

510 

1010 

1968 

3682 

6222 

8072 

5126 

492 

0 

16 

32 

64 

511 

1018 

2013 

3910 

7268 

12112 

15239 

9044 

750 

0 

B 

32 

64 

128 

256 

1022 

2034 

|||||[QQg 

7776 

14362 

23608 

28824 

16028 

1164 

0 

Qj 

256 

512 

1023 

2042 

4061 

8006 

B 

B 

128 

256 

512 

1024 

2046 

1  4082 

8112 

15968 

30728 

B 

B 

512 

1024 

8157 

16198 

31836 

61074 

110634 

174507 

4305 

B 

512 

16304 

32352 

63490 

121422 

218542 

340034 

368356 

158598 

6710 

B 

26 

1024 

llglg 

32582 

64604 

126592 

241352 

431607 

662419 

696663 

281830 

0 

B 

16370 

32688 

65120 

129024 

252442 

852524 

1290716 

1318138 

16264 

0 

Bl 

4096 

130140 

257658 

B 

29 

8192 

260096 

514568 

B 

B 

16383 

519800 

1027594 

2001434 

3768843 

6568697 

9547486 

8933851 

2837467 

61967 

0 

B 

32766 

1038850 

2052174 

3990960 

7492056 

1.30E+07 

1.68E+07 

1.68E-K)7 

5064334 

97108 

0 

Table  21,  «  =  31  (part  I) 
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